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Electromagnetic Oscillations in a Driven Nonlinear Resonator: 
A New Description of Complex Nonlinear Dynamics 

E. Yu. Petrov and A. V. KudrirE 

Department of Radiophysics, University of Nizhny Novgorod, 

23 Gagarin Ave., Nizhny Novgorod 603950, Russia 

Many intriguing properties of driven nonlinear resonators, including the appearance of chaos, are 
very important for understanding the universal features of nonlinear dynamical systems and can 
have great practical significance. We consider a cylindrical cavity resonator driven by an alternating 
voltage and filled with a nonlinear nondispersive medium. It is assumed that the medium lacks a 
center of inversion and the dependence of the electric displacement on the electric field can be ap- 
proximated by an exponential function. We show that the Maxwell equations are integrated exactly 
in this case and the field components in the cavity are represented in terms of implicit functions 
of special form. The driven electromagnetic oscillations in the cavity are found to display very in- 
teresting temporal behavior and their Fourier spectra contain singular continuous components. To 
the best of our knowledge, this is the first demonstration of the existence of a singular continuous 
(fractal) spectrum in an exactly integrable system. 

PACS numbers: 05.45.-a 



Nonlinear resonators are simple and convenient mod- 
els of physical systems and have proven to be a valuable 
means to investigate the universal features of nonlinear 
dynamics. In view of this, such resonators have been 
the subject of intense theoretical and experimental stud- 
ies in the past decades [l|-|j]. A variety of the existing 
electronic devices and materials with nonlinear electro- 
magnetic properties makes it possible to create electrical 
resonators with different types of nonlinearity. From a 
practical viewpoint, it is important that when electro- 
magnetic oscillations in a bounded volume are excited at 
frequencies close to resonant ones, oscillation amplitudes 
that are sufficient for manifestation of various nonlinear 
phenomena can be reached even if a driving source is 
comparatively weak. The well-known example of such 
phenomena is generation of harmonics of a drive fre- 
quency, which enables realization of efficient frequency 
multipliers on the basis of high-Q resonators. In addition, 
fairly complex, e.g., chaotic, oscillations can be excited 
in nonlinear resonators [l|-|3|. Studying chaotic regimes 
of resonators is of great interest for developing control- 
lable sources of noise-like signals. In more recent works, 
a new type of complex nonlinear dynamics, intermediate 
between almost periodic and random, has been discov- 
ered [5|-|7|. Such dynamics, which is associated with a 
singular continuous spectrum, appears typically in driven 
nonlinear systems and has received much attention [8|- 
[14] , both for fundamental reasons and because of its in- 
terdisciplinary relevance [T'l . 

A complete description of the complex dynamics of dis- 
tributed nonlinear systems is fairly difficult to achieve. 
This is explained by an infinite number of degrees of free- 
dom and the presence of several controlling parameters 
in such systems. Because of this, most theoretical works 
on the subject discuss nonlinear resonators as lumped 
systems or merely as a collection of coupled oscillators 



or modes. Within the framework of such a simplified 
approach, the problem of oscillations in a nonlinear res- 
onator is reduced to solving a system of ODEs. Although 
such an approach is justified in many cases, it is clear that 
electromagnetic systems should generally be described by 
the Maxwell equations. 

In this work, the problem of a nonlinear electrical res- 
onator is considered using a full set of the Maxwell equa- 
tions. In what follows, we apply the method for con- 
structing exact axisymmetric solutions of the Maxwell 
equations in a nonlinear nondispersive medium, which 
has been developed in our recent works 151, ll6| , to the 
driven oscillations in a bounded volume. 

Consider electromagnetic fields in a cylindrical cavity 
of radius a and height L. We assume that the z axis of a 
cylindrical coordinate system (r, (f), z) is aligned with the 
cavity axis and limit ourselves to consideration of axisym- 
metric field oscillations in which only the E^ and H^ com- 
ponents are nonzero. We will also assume that the cavity 
is filled with a nonlinear nondispersive medium in which 
the longitudinal component of the electric displacement 
can be represented as D^ = Dq + a^^eoffi [e:xjp{aEz) — 1], 
where Dq, ei, and a are certain constants. The pos- 
sibility of using such a model of nonlinearity for media 
lacking a center of inversion is discussed in detail in [l5| . 
As is shown in [15| . this model, with appropriately chosen 
Do, £i, and a, correctly describes dielectric properties of 
such media in the case of moderately small electric fields. 
Then the Maxwell equations read 



drH 



'H = eiE)dtE, drE = fiodtH, (1) 



where E = Ez{r,t), H = H^{r, t), and 

e{E) = dDjdE = eoSi exp{aE). (2) 

The exact solution of system (1) can be written in im- 



plicit form as [15J 

E = A-^E (p e"^/^ T + apH/2 

H = e"^/2^-iH (p e'^^'^T + apH/2 



(3) 



Hereafter, A is a constant amplitude factor related to the 
field source, E = E/A, H = ZQH/{Ae\''^), p = r/a, t = 
t{eoeip,o)^^^^/a,, and a = aA, where Zo = (/xo/eo)"^^^- 
The functions £ and H describe the electromagnetic field 
in a linear medium and satisfy the equations 



d^S 



-'dp£ 



dl£, 



(4) 



and dpE 



drU. 



Let us state the following initial and boundary condi- 
tions for linear wave equation (4): 



(5) 



£:(/9,o) = o, arf(p,o) = o, o<p<i, 

£(1,t) =yl(sinfiiT + 0.5 sinfiar), < r < oo, (6) 

where r2i_2 are normalized constant frequencies such that 
f^i,2 = Ai,2a(eo£iAio)^^^ (i-e., ^i,2T = Xia^)- The fre- 
quencies f2i and O2 are related as fii = <7^2^ where 
(7 — (Vs — l)/2 is the golden mean. The boundary value 
problem defined by Eqs. (H])-® describes the driven elec- 
tromagnetic oscillations in a cylindrical cavity specified 
by the relations p — r /a < 1 and < 2: < L, which 
is filled with a linear medium having the permittivity 
e ~ eo^i = const {a = 0). The linear oscillations are 
thus driven by an electric field ^ which can be pro- 
duced by two coaxial metal rings of radius a that are 
separated by distance L if an almost periodic voltage 
V = £{\,t)L — AL{sv[v\it -f 0.5 sinA2t) is applied be- 
tween them. The solution to the linear boundary value 
problem specified by Eqs. (HJ-® can be found in a stan- 
dard way [1^. As a result, the functions E and % are 
written as 



£:(p,r) = ^B, 



2 






Join,) 



sin fljT- 



cos flj T - 



n=\ 



CnJoiKnp) sin KriT, 



y^ C„Ji(k„p)cosk„t, 

n=l 

(7) 



where B, = A/j, C„ = 2Aj:^^^, n,[j{n]-Kl)J,{Kn)]-\ 
Jm is the Bessel function of the first kind of order m, k„ 
is the nth root of the equation Jo(k) = 0, and fli 2 7^ ««■ 
We denote eigenfrequencies of the E'ono modes [l8| of a 
linear resonator as w„. Hence, aj„ = K„(eo/^o)~^ a~^ 
and K„T — ujnt. 

Substituting functions ^ into formulas ([3]) , we obtain 
an exact solution to system ([T]) in implicit form. Thus, 
the field components E and i7 in a cylindrical cavity 
filled with a nonlinear medium for which £{E) is written 
in the form of Eq. ^ are found from the solution of 



a set of transcendental equations ((S)) in which £ and T-L, 
defined by relationships ^ , are almost periodic functions 
of r. In the nonlinear case, the fields E and if, which are 
determined by Eqs. (|3]) and ([7]), satisfy the same initial 
conditions as in Eq. (5): 

£;(p,0)==0, drE{p,Q)^Q, 0<p<l. (8) 

However, the electric field oscillations on the side wall of 
the nonlinear resonator (p = 1) in the found solution do 
not obey Eq. (6). Putting p := 1 in formulas (3), we have 

iB/2 



E = A-^S (e" 



aH/2 



H 



E/2 



&H/2 



(9) 



The dependence E{\,t) determined by relationships (|9]) 
can be regarded as a drive signal in the boundary value 
problem given by Eqs. ([1]), (jH)), and ^ for a nonlinear 
(a 7^ 0) resonator. 

Thus, formulas ([3|), with E and T-L given by relation- 
ships ([7]) , yield an exact solution of the nonlinear bound- 
ary value problem for system ^ under conditions ([5]) 
and (|n]), and describe driven electromagnetic oscillations 
in a cylindrical cavity resonator filled with a nonlinear 
medium. A typical example of such a medium can be 
a ferroelectric crystal. Note that ferroelectric resonators 
are known to be of great interest for many promising 
applications [J]. 

Observe that the oscillations on the axis p = of the 
nonlinear resonator in the obtained exact solution coin- 
cide with the oscillations for p = in the "seeding" linear 
problem. It follows from Eqs. ([3|) and ([7]) that for p = 0, 
-B(0, r) = f (0, r) = Y.]=i Bj sinrj^r + Yln=i ^n sinK„T 
and i?(0, r) = H(0, r) = 0. Thus, the electric field oscil- 
lations on the axis p = are described by an almost pe- 
riodic function of t and have the discrete spectrum [19J 1 . 

For p 7^ 0, the exact solution expressed in terms of 
implicit functions is very complicated and can be studied 
only numerically. It turns out that for p 7^ 0, the field 
oscillations described by this solution may have singular 
continuous (fractal) spectrum. 

To confirm the above assertion, we turn to results of 
calculations of the quantities E and H determined by 
Eqs. ([3|) and ([7]). In what follows, the main attention 
will be focused on analyzing the obtained solutions in 
the case where the drive frequency A2 relates to the 
fundamental eigenfrequency uii of a linear resonator as 
^2/^^! = £^2/^1 = a (note that the identity VL1 + Q.2 = ni 
is valid in this case). In our calculations, we retain 100 
terms of the series over n in formulas ([7]). It should be 
noted that the employed theoretical model of a nondis- 
persive medium does not allow one to indefinitely in- 
crease the nonlinearity parameter a in solutions (3). For 
large absolute values of a, implicit functions E{p, r) and 
H{p,t) determined by Eqs. ([3]) and (O become ambigu- 
ous, and solutions ([3]) obtained without allowance for dis- 
persion will be inapplicable [15]. For fixed a and r2i^2, 



the ambiguity points appear first in the time dependences 
E{1,t) and H{1,t) (for p = 1) since the nonlinear effects 



are accumulating with increasing p 15|. Therefore, as a 
first step in practical calculations, one should study the 
functions E{1,t) and H{1,t). If these functions are un- 
ambiguous and continuous, then E and H possess the 
same properties for < p < 1. In all the computations, 
we use the maximum possible value a = 0.32 for chosen 

Now consider the field oscillations on the wall p — 1 oi 
the nonlinear resonator. The dependences E{1,t) and 
H(1,t) determined by relationships ([9]) are shown in 
Fig. 1(a) by the red and blue lines, respectively. For 
comparison, the solid and dashed black lines in Fig.l (a) 
show the functions £ and H, respectively. It is seen in 
Fig. 1 that the functions E and H demonstrate fairly 
complex behavior and essentially differ from the corre- 
sponding quantities £/A and H/A in the linear regime 
{a = 0) by the presence of small amplitude spikes. Fig- 
ures 1(b) and 1(c) show the frequency spectra of E{1,t) 
and H{1,t) (2^^ points were taken for the FFT with a 
sampling rate of 0.05). The FFT plots exhibit numerous 
spectral components that can be attributed to the har- 
monics of rii.2 and K„, as well as various combination 
frequencies. 

To reveal nontrivial spectral properties of the oscil- 
lations considered, we use the singular-continuous spec- 
trum analysis |6l.l9l4lll.ll3|. We define the partial Fourier 
sums 



SEin,T)= J2e„ 



T 
m=l 



(10) 

where {Em} and {H„i} are the time series of the variables 
E and H: Em = E{Tm) and Hm — H{Tm)- We take n — 
100 and Tm+i — Tm — 0.02, although all the forthcoming 
results are valid when changing the initial point and the 
sampling rate. The Fourier transforms scale with T as 



|5B(r!,T)|2^T^, \SH{n,T)\'' ^T\ 



(11) 



where 5 = B(^) and 7 — 7(ri) are scahng expo- 
nents [6|, y, lul . The evolution of Se and Sh with T can 
be represented by paths in the complex planes (ReS'^;, 
\t[iSe) and (ReS'ij, ImS/f), respectively. It is known [9] 
that for /3 = 7 = 2, the frequency SI belongs to the 
countable set of discrete spectral components of an al- 
most periodic oscillation and there exist persistent mo- 
tions (drifts) of Se and Sh in the corresponding com- 
plex planes. A singular-continuous spectral component 
appears if (i) (3 ^ 1,2 and/or 7^ 1, 2 and (ii) the path 
in the complex plane is fractal [y, I9l4lll Il3| . A singular 
continuous spectrum is known to be a Cantor set [71, [ij] . 
We have found that at some frequencies, the spectrum 
has the scaling exponents /3 = 7 = 2. Figures 2(a) and 
2(b) show \Se\'^ and \Sh\'^ as functions of logj^pT for p = 1 
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FIG. 1. (color online) (a) Oscillograms of the electric (E) and 
magnetic (H) fields on the wall p = 1 of a nonlinear resonator 
(red and blue lines, respectively), calculated by formulas (9), 
and the corresponding quantities £/A and H/A in the linear 
regime (solid and dashed black lines, respectively). The fast 
Fourier transforms of (b) E and (c) H . 



at one of such frequencies, namely, fl — k^ = 8.65 .... 
The corresponding paths in the complex planes are pre- 
sented in Figs. 2(c) and 2(d). Thus, in this case, we deal 
with a discrete component of the spectrum. However, at 
the combination frequency ^1,7,3^4 = 32.96 . . . (hereafter, 
^i,m,n = l^i+mQ2 + nK2), wc havc /3 ~ 1.7 and 7 ~ 1.62 
[see Figs. 2(a) and 2(b)], and the behavior of the depen- 
dences IS'sp and l-Sj? P in this case is typical of a singu- 
lar continuous component [9|, IllL |l3| . The corresponding 
paths in Figs. 2(e) and 2(f) exhibit fractal structures. 
These results strongly suggest that the considered spec- 
trum of electromagnetic oscillations is not purely discrete 
and contains singular continuous components. For p = 1, 
we have also found such components at, e.g., the combi- 
nation frequencies 03^3^4 (/3 « 1.57 and 7 = 0), ^23,44 
(P « 1.6 and 7 = 0), '^^6,3,4 (/3 ~ 1-53 and 7 « 1.4), 
f^3,-i,4 (/3 w 1-1 and 7 « 1.55), ^^6,5,4 (/? » 1-52 and 
7 w 1.8), and ri5,6,4 (/3 ~ 1-4 and 7 ~ 1.2). For many 
frequencies, a power-law growth of the spectrum is ob- 
served with the exponents /3 and 7 which differ from 2 
only slightly. 

Wc now pass to consideration of some spectral features 
of oscillations inside a nonlinear resonator for p = 0.5. 
Here, the components of the singular continuous spec- 
trum appear at higher frequencies than for p = 1. The 
values of /3 and 7 for p = 0.5 turn out to be smaller 
than for p = 1 at the same frequency. For example, at 
il = rig. 5, 4 — 35.02 . . ., we have /3 w 1.52 and 7 w 1.8 
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FIG. 3. (color online) Singular-continuous spectrum analysis 
of the time series {Em\ and {ilm\ for p — 0.5. (a) \Sb^ 
and (b) \Snf' as functions of logj^gT aX Q. — f26,5,4. The 
corresponding paths of (c) Sb and (d) Sn in the complex 
planes. 
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FIG. 2. (color online) Singular-continuous spectrum analysis 
of the time series {-&,„} and {-ff„i} for p = 1. (a) IS'sp and 
(b) \Sh^ as functions of log^T at f2 = K3 and Q. — f27,3,4. 
The paths of (c) Se and (d) Sh in the complex planes (Re Se , 
IhiSe) and (RbSh, ImSH), respectively, at 57 = K3. The 
paths of (e) Se and (f) Sh in the same planes at f2 = fir, 3, 4. 



for p = 1, and ^ ~ 1.05 and 7 « 1.05 for p = 0.5 
[see Figs. 3(a) and 3(b)]. The corresponding curves in 
the complex planes (RbSe, IthSe) and (ReSn, ImS*//), 
which are presented in Figs. 3(c) and 3(d), display fractal 
behavior. 

Thus, our analysis shows that the Fourier spectrum of 
the electromagnetic oscillations in the cavity is a mixture 
of discrete and singular continuous components. Simi- 
lar phenomena have been reported in the literature and 
discussed in, e.g., [9|, [ij] as applied to the dynamics de- 
scribed by forced maps and symbolic sequences. In our 
case, calculations of the FFT and the autocorrelation 
function (ACF) do not allow one to adequately investi- 
gate the spectral properties of the considered oscillations. 
Because of the weakness of the contribution from the sin- 
gular continuous part, the FFT and the ACF are similar 
to those of almost periodic motion. Scaling is a powerful 
method to detect the presence of a singular continuous 
spectrum. However, it is not clear how this component 
can be separated from the discrete spectrum. 

The existence of regimes with singular continuous 
(fractal) spectra in dissipative dynamical systems de- 
scribed by discrete maps or ODEs is well known |8l4ll|. 
Such regimes corresponding to strange nonchaotic attrac- 
tors are realized on sets of positive measure in the param- 
eter spaces of dissipative dynamical systems and are typ- 



ical of the intermediate region between almost periodic 
and random motions. Our study demonstrates that the 
nonlinear dynamics with a singular continuous spectrum 
can occur in an exactly integrable distributed nondissi- 
pative system. We have found that the implicit functions 
given by Eqs. ([3]) and (O, which are exact solutions of 
system ([T]) , are not almost periodic in r and their Fourier 
spectra contain singular continuous components. Study- 
ing these functions is of great interest for physics and 
mathematics. We have shown that such functions can be 
finite-amplitude single-valued solutions of the Maxwell 
equations and, hence, describe actually existing electro- 
magnetic oscillations. Thus, Eqs. (jS)) and ^ provide a 
new description of complex nonlinear dynamics. 
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